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ABSTRACT 

In this work, we examine the different properties of galactic satellites in hydro- 
dynamical and pure dark matter simulations. We use three pairs of simulations (col- 
lisional and collision-less) starting from identical initial conditions. We concentrate 
our analysis on pairs of satellites in the hydro and Nbody runs that form from the 
same Lagrangian region. We look at the radial positions, mass loss as a function of 
time and orbital parameters of these "twin" satellites. We confirm an overall higher 
radial density of satellites in the hydrodynamical runs, but find that trends in the 
mass loss and radial position of these satellites in the inner and outer region of the 
parent halo differ from the pure dark matter case. In the outskirts of the halo (« 70% 
of the virial radius) satellites experience a stronger mass loss and higher dynamical 
friction in pure dark matter runs. The situation is reversed in the central region of 
the halo, where hydrodynamical satellites have smaller apocenter distances and suffer 
higher mass stripping. We partially ascribe this bimodal behaviour to the delayed in- 
fall time for hydro satellites, which on average cross the virial radius of the parent halo 
0.7 Gyrs after their dark matter twins. Finally, we briefly discuss the implications of 
the different set of satellite orbital parameters and mass loss rates in hydrodynamical 
simulations within the context of thin discs heating and destruction. 



Key words: galaxies: haloes - cosmology:theory, dark matter, gravitation - 
numerical, N-body simulation 



methods: 



1 INTRODUCTION 

In the current paradigm of structure formation, large ob- 
jects, such as galaxies or clusters, are believed to form hier- 
archically, through a 'bottom-up ' (|White fc Reesll 19781 ) pro- 
cess of merging. About a decade ago, N-body simulations 
attained sufficient dynamic range to reveal that, in Cold 
Dark Matter (CDM) models, all haloes should contain a 
large number of embedded subhaloes that su rvive the col- 
lapse and virialization o f the parent structure l|Klypin et al.l 
ll999l : lMoore et alJll999h . 

The properties of subhaloes on different scales has been 
the subject of many recent studies that ha ve pushed the 



resolution of dissipationles s simulations (e.g. ISpringel et al 
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De Lucia et al.ll2004l: iKravtsov et al.ll2004 



Gao et al 



2 0051: IDiemand et al.ll2007l ; IZentner et al 
20081 ) 



The kinetic properties of sub- 
haloes are now well understood - they make up a fraction 
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of between 5 and 10% of the mass of virialized haloes, on 
scales relevant to observational cosmology. 

Most of these previous studies used dissipationless cos- 
mological simulations; although non-baryonic dark matter 
exceeds ba ryonic matter by a fa ctor of Sldm/Qb — 6 on av- 
erage (e.g. iKomatsu et~aiTl2009l ). the gravitational field in 
the central region of galaxies is dominated by stars and gas. 
The cooling baryons increase the density in the central halo 
region mainly because of the extra mass associated with the 
inflow, but also because of t he adiabatic contra ction of the 
total mass distribution (e.g. iGnedin et al"1l2004l ). Since this 
process is active for both the host halo and its subhaloes, it 
might be expected that subhaloes formed within hydrody- 
namical simulations (including gas and stars) will experience 
a different tidal force field and will themselves be more ro- 
bust to tidal effects. 

Recently, a number of authors have examined the 
impact of baryonic physics (gas cooling, star formation 
and feedback) on both the central object and the satel- 
lite populat i on in galaxy and clust e r-size d haloes (e.g. 
I Baffin et al. | 120051 : iNagai fc Kravtsovl 120051 : iMaccio et all 
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2006: 



Weinberg et alJ 



20081: iRomano-Dfaz et aHl2009l . [201ol ; 
Libeskind et al.l I201C ; ISomme r-Larscn fc Limousi 3 l20ld ; 
Duffy et al.ll2010t ). IMaccio et all l|2006h simulated a Galac- 



tic mass halo twice - once considering pure DM and once 
including baryons modeled with smoothed particles hydro- 
dynamics, stopping the gas cooling at z = 1.5. They found 
that the hydro run produced an overabundance of subhaloes 
in the inner regions of the halo as well as an increase by a 
factor of 2 in the absolute number of subhaloes with respect 
to the DM run. 

The issue of the distribution and properties of galac- 
tic satellites in hydro and DM simulation s has been 
more recently revisited bv iLibeskind et al] (fe oiO'l and 
iRomano-Di'az et al l (|201Gi ). ILibeskind et al l (|2010|7 found re- 
sults very similar to the work of Maccio et al.l ( 20061 ). with 
subhaloes in the hydro simulation being more radially con- 
centrated than their dark matter counterparts. They ascribe 
this effect to the higher central density of subhaloes in hy- 
dro simulation (due to the collapse of baryons into stars in 
the central region) that makes them more resilient to tidal 
forces. The increased mass in hydrodynamic subhaloes with 
respect to dark matter ones causes dynamical friction to be 
more effective, dragging the subhalo towards the centre of 
the host. 

The overall properties o f the s atellite population in the 
study of Romano-Diaz et al. 1 20icil) are al s o con sistent with 
IMaccio et al.1 (|2006l ) and ILibeskind et al l (|201Cl ). But, per- 
haps counter-intuitively, they find that satellites in the hy- 
dro run are depleted at a faster rate than the pure DM one 
within the cen tral 30 kpc of the p r ime h alo. According to 
their analysis, IRomano-Di'az et al.l (|20101 ) suggest that al- 
though the baryons provide a substantial glue to the sub- 
haloes, the main halo exhibits the same trend. This would 
assure a more efficient tidal disruption of the hydro subhalo 
population in the inner region of the halo (~ 0.1 of the virial 
radius) . 

Studies concerning the different results of pure dark 
matter and hydrodynamical simulations are especially com- 
pelling within the context of satellite effects on disk sta- 
bility. Several theoretical and numerical studies have been 
devoted to quantifying the resilience of galactic disks to 
infalling satellites (e.g. iToth fc Ostrikerlll992l; Quinnetal 



Table 1. Galaxies parameters 



Galaxy 


Mass 


Rvir 


Nsat 




(10 11 h- 1 M e ) 


(kpc/h) 




GO (DM/Hydro) 


7.8/7.4 


191/188 


63/77 


Gl (DM/Hydro) 


9.2/8.9 


201/199 


70/87 


G2 (DM/Hydro) 


10.0/9.3 


207/202 


102/110 



2008 



1993 1; I Velazquez fc Whitell 19991: 



Font et al.ll2001l : |Read et al. follow ing lHaardt fc Madaul 
Moster et al.1 120101 ). Re- 



radial position, and orbital parameters (peri and apo-centre 
distances). We start hydro and pure dark matter simulations 
from the same initial conditions and we focus our analysis on 
"twins" satellites, i.e. (sub)structures that are formed from 
the same Lagrangian region for the initial conditions, which 
should therefore share the same formation history in both 
simulation types. 

The remainder of the paper is organized as follows: in 
Section [5] we describe our simulations, and provide a brief 
summary of the numerical codes we use, including the tech- 
nique employed to match satellites in the different runs. In 
Section [3] we present our main results, focusing on several 
satellite properties like radial position, orbital parameters, 
mass loss. Finally in Section [4] we summarize and discuss 
our results. 



2 NUMERICAL SIMULATIONS 

The hydro- simulations were performed with GASO- 
line, a multi-stepp ing, parallel TreeSPH iV-body code 
|Wadslev et al.l [20041 ). We include radiative and Compton 
cooling for a primordial mixture of hydrogen and helium. 
The star formation a lgorithm is based on a Jeans instability 
criteria (|Kata Il992j ). but simplified so that gas particles 
satisfying constant density and temperature thresholds in 
convergent flows spawn star parti cles at a rate propor tional 
to the local dynamical time (see IStinson et al.l [2006) . The 
star formation efficiency was set to 0.0 5 based on simu- 
lations of the Milky Way satisfying the iKennicuttl (| 1 9981 ") 
Schmidt Law. The code also i nclud es supernova feedback in 
the manner of Stinson et al.l 



Villalobos fc Helmil |2008[ 

Kazantzidis et al.1 (|2009T ) performed a very detailed 



2006), and a UV background 



1996) (see l|Governato et al.l 



cently 

study of the dynamical response of thin galactic disks to 
bombardment by cold dark matter substructure. They used 
pure Nbody simulations of the formation of a Milky Way- like 
dark matter halo to derive the properties of substructures 
and subsequently as initial conditions in subsequent high 
resolution satellites-disk merger simulations. Clearly, under- 
standing i/and how different possible orbital parameters and 
mass loss rates expected for satellites in hydro runs modify 
results previously obtained using pure Nbody simulations is 
of importance. 

In this work we revisit the issue of the effects of baryonic 
physics on the satellite population in galactic-s ize dark mat- 
ter ha loes. We improve the original study bv IMaccio et al.l 
(2006) in several aspects, namely with better parametriza- 
tion of the baryonic physics, a full hydrodynamical approach 
down to redshift zero and a more extensive analysis of the 
satellite properties, including the time evolution of mass loss, 



l2007f) for a more detailed description of the code). 

Three candidate haloes with masses similar to the mass 
of our Galaxy (M ~ 1O 12 M ) were selected from an exist- 
ing low resolution dark matter simulation (300 3 particles 
within 90 Mpc) and subsequently re-simulated at higher 
resolution. These high resolution runs are 8 3 times more 
resolved in mass than the initial set and we included a 
gaseous component within the entire high resolution re- 
gion. Masses of the dark matter and gaseous particles are 
m d = 1.17 x 10 6 /i -1 M and m 9 = 2.3 x W 5 ^ 1 M , respec- 
tively. The dark matter has a spline gravitational softening 
length of 500 h -1 pc and each component (dark and gas) 
consists of about 10 6 particles in the high resolution region. 
A list of galaxies properties can be found in Table [1] A more 
detailed description of these hydro simulations (with partic- 
ular attention to the properties of the central galaxy) will 
appear in a forthcoming paper (Hernandez et al. in prep). 

To run the dark matter onl y counter-par ts of G0-G2 we 
use the Nbody code pkdgrav ((stadcl 2001). This code is 
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Figure 1. Radial density profile for G0-G2 in the hydro (red 
lines) and DM (blue lines) simulations. 



intimately related with GASOLINE which is its hydrodynam- 
ical extension. This implies that DM particles are treated in 
exactly the same way in the two codes (e.g. using the same 
numerical algorithms), allowing a straight, direct compari- 
son of the results. The initial conditions for the DM sim- 
ulations are identical to those used for the hydro runs; we 
simply transform all gas particles into dark matter particles. 
This insures us both that the normalization and the phases 
of the initial density perturbations are identical in the two 
runs. These Nbody simulations are the same presented in 
iMaccio et all (|2010l ). 

Figure [T] shows the density radial density profiles of the 
G0-G2 galaxies in the two runs: hydro (red lines) and DM 
(blue lines). As expected the profiles only diverge in the 
central regions due to the presence of a baryonic core in the 
hydro runs. 

2.1 Halo finder and merger tree construction 

To identify subhaloes in our simulation we use the 
MPI+OpenMP hybrid AHF halo finder, available freely at 
http://www.popia.ft .uam.es /AMIGA ) and described in de- 
tail in Knollmann & Kncbc (2009). AHF identifies local over- 
densities in an adaptively smoothed density field as prospec- 
tive halo centers. The local potential minima are computed 
for each of these density peaks and the gravitationally bound 
particles are determined. Only peaks with at least 50 bound 
particles are considered to be haloes and retained for further 
analysis. As subhaloes are embedded within their respective 
host halo, their own density profile usually shows a charac- 
teristic upturn at a radius rt < r v i r , where r v i r is the actual 
(virial) radius of the satellites, were they found in isolation. 
We use this "truncation radius" rt as the outer edge of the 
subhaloes and calculate all (sub-)halo properties (i.e. mass) 
using only the gravitationally bound particles inside rt. 

As second step we build the merger trees for our galaxies 
(both hydro and DM), in order to establish the dynamical 



history of each infalling satellite across cosmic time. For the 
purpose of constructing an accurate tracking of each simu- 
lated (sub)halo, we analyse 40 simulation outputs between 
z — 2.5 and z — 0. We start from the halofinder results at 
z — 2.5 and follow these (sub)haloes through cosmic time, 
adding to the tracking procedure at each snapshot all new 
haloes that cross our mass threshold. For halo tracking, only 
dark matter particles are used. We consider two criteria to 
decide if, when comparing two consecutive snapshots, halo 
1 at one output time is the "progenitor" of halo 2 at the 
subsequent time step: i) more than 50% of the particles in 
halo 1 that end up in any halo at time step 2 end up in halo 
2; ii) more than 50% of the particles in halo 2 come from 
halo 1. 

2.2 "Twin" matching and orbital parameters 

Since we aim to compare the hydrodynamical and pure DM 
simulations by studying satellite "twins", i.e. satellites that 
formed in both simulations from the same perturbations at 
the initial time, we perform the following steps to identify 
these satellite pairs: 

(i) We use the AHF halo finder to identify a given satellite 
in the DM simulation at the redshift z of interest. 

(ii) We track the position of the satellite DM particles 
back in time to the initial condition file, in order to deter- 
mine its original Lagrangian region. 

(iii) We look for DM particles within the set of hydro 
initial conditions in the same Lagrangian region and map 
them forward in time to the same redshift z. 

(iv) Using the satellite catalogue for the hydro simulation 
we check in which satellite(s) these dm particles end up. 
We then use the same criteria adopted in the merger tree 
construction to determine the correct "twin" pair (see Fig. 

ED. 

We have explicitly tested that the exact same results for 
the satellite mapping follow starting either from the DM or 
the hydro simulation. Overall we find more than 200 satel- 
lites pairs and 76 of them survived down to z = in both 
simulations. We will refer to this last subset of satellites as 
the twin population at redshift zero. 

Finally we compute the orbital parameters (apo-center 
and peri-center distances) for each pair of twin satellites. 
In order to do that we integrate the satellite orbit in an 
effective static potential paramet rized using the powe r-law 
logarithmic slope (PoLLS) model (|Cardone et alj|2005l ): 



polls / \ / 2 
P (r) = p_ 2 exp 

7 



R- 



- 1 



(1) 



where p_2 is the density at the R-2 distance, which is where 
the profile slope is equal to —2 (i.e. the singular isother- 
mal sphere profile) and 7 is related the central slope of the 
profile. All these parameters are fitted to match the actual 
potential in either the DM or hydro simulation. In prin- 
ciple since we have constructed the merger trees it would 
have been possible to compute the satellite orbital parame- 
ters directly from the simulation. On the other hand, given 
the time sampling of our snapshots we could end up with 
over/under estimating the peri/apo-centre distances. So we 
decided to use a frozen potential in order to have a more 
reliable orbital parameter estimate. 
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Figure 2. Illustration of satellite mapping scheme. 



3 RESULTS 

In this section we study how the presence of baryons affects 
the properties of satellites in our simulations. We include in 
our analysis only substructures with at least Mtot > 1O 8 M0. 
We combine together the results of GO, Gl and G2 (which 
give similar results when analyzed individually) and mainly 
focus on "twin" haloes. This gives us a sample of more than 
70 satellite couples at z — 0. In the following we refer to 
the dissipationless simulations as pure DM or Nbody, while 
we use the therm "hydro" for dissipational simulation that 
include gas and stars. 

Figure [3] shows the cumulative radial distribution (for 
all satellites) in our three galaxie s in the two d if ferent sim- 
ulation runs. As first noticed by iMaccio et al.l (120061) and 
later confirmed by other studies (e.g. iLibeskind et al.ll2010l : 
iRomano-Diaz et al.ll2010T ) the radial distribution of satellites 
is clearly more concentrated in hydro simulations (solid red 
line) than in DM ones (dashed blue line). This trend is par- 
tially confirmed when we restrict our analysis to twin satel- 
lites as shown in Figure [4j at least in the inner region. We 
should point out that "twin" satellites are a biased sub sam- 
ple of the whole satellite population, since, by construction, 
we restrict ourselves to those satellites that survived until 
z = in both simulations. This biased selection explains 
the difference in the radial distribution profiles of the whole 
population and of the twin sample in the hydro and DM 
runs (especially in the external regions). 

In order to understand the physical drivers behind the 
results of Fig. [3] and Fig. 0] we exam individual mass accre- 
tion and dynamical histories of all twin couples from z = 1.2 
to the present time. This is the redshift at which the latest 
major merger happen (namely for the Gf halo), after z = 1.2 
our galaxies have a quiet dynamical history that allows to 
study in details the satellites properties and orbits Figure [5] 
shows the evolution with redshift of the distance from the 
halo centre (upper panel) and the total mass (lower panel) 
of an hydro subhalo (solid line) and its dm twin (dashed 
curve). The two satellites enter the virial radius of the main 
object (R = 1) at the same time and have similar orbits, 
ending at z — at roughly the same distance from the halo 
centre. In this particular case the mass loss is higher in the 
hydro case, where the satellite is able to retain only 30% of 
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Figure 3. Number density profile of ALL satellites at z = 0. 
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Figure 4. Number density profile of TWIN satellites at z = 



its mass (50% in t he pure DM run) . This figure is similar to 
Figure 8 and 9 of ILibeskind et all ||2010I l even though the 
role of DM and hydro are reversed. 

The behaviour presented in Fig.[5]is rare, however. The 
majority of satellites have more complicated dynamical his- 
tories as shown in Fig. [6] There, it is immediately noticeable 
that the twins in this case do not enter the virial radius of the 
parent halo at the same time. The infall time is z in f a u = 0.93 
in the hydro case and z in f a ii = 1.25 in the pure DM one. The 
orbits are moreover very different: in the hydro simulation 
the satellite has a very large apo-center that brings it out- 
side the virial radius at z ~ 0.5 before being re-accreted; this 
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Figure 5. Evolution of radial distance (upper panel) and mass (lower 
panel) over time for a satellite that lost more mass in the hydrodynamical 
simulation than in the pure DM one. 
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Figure 6. Same as figure^ but for a satellite that lost less mass in the 
hydrodynamical run than in the pure DM one. 



does not happen in the pure DM run. As a consequence, the 
orbits are practically uncorrelated, and satellites experience 
different mass loss rates. 

An interesting feature revealed by Fig.|6]is the difference 
in the time of accretion between the hydro and DM runs. 
This is not a peculiar behaviour of the selected satellite, but 
it is quite a general trend, as shown in Fig. [7] For each twin 
couple we compute the infall time (defined as the time at 
which the satellite first crosses the virial radius of the main 
object) in the hydro and pure DM runs and then we plot 
the difference in these two times as a function of the infall 
redshift in the DM case. 

Most of the hydro satellites enter the virialized region 
later than their pure DM counterparts, with an average de- 
lay of 0.7 Gyr^j. A possible explanation for this delayed 

1 Some haloes show a very large At. These are twins were one of 
the two subhaloes went in and out the viri al radius before be ing 
accreted. The so-called back splash haloes l|Knebe et alj|2010r i 



Figure 7. Difference in infall time (e.g. crossing of the virial radius) for 
twins. ti n f a u is the time measured from the Big-Bang. 



accretion might be time evolution in the main halo virial 
radius. Fig. [9] shows R V i r (z) for the hydro and pure DM 
simulation of the G2 galaxy. (We observe a similar effect for 
GO and Gl also.) In both runs the virial radius (in comoving 
coordinates) grows until z w 0.8 and then decreases toward 
z = 0. Although the DM run shows a larger value of R V ir 
after z = 1, this difference is very small (about 5%) and can- 
not explain the results of Fig. [7] This also confirm by figure 
[5] where we plot the difference of physical distance from the 
center of the halo for twin couples at the moment the hy- 
dro satellite crossed the virial radius. The plot confirms that 
hydro satellites are substantially farther away compared to 
their DM counterparts at the time of accretion. 

An alternative, possible, explanation arises from the 
point of view of the initial conditions. The hydro and pure 
DM runs differ only in that a fraction ttb/Qdm of the total 
mass in the hydro simulation is represented by gas particles, 
while these gas particles are converted to dark matter par- 
ticles for the Nbody run. The lower panel of Fig. [9] shows 
the mass within a fixed radius of 270 kpc due to baryonic 
particles (stars + gas) in the hydro run and the same DM- 
converted particles in the pure DM run. After z = 1 there, 
the amount of "converted particles" inside a sphere of 270 
kpc exceeds the baryonic mass accumulated within the same 
region. One possible explanation for this effect is that the 
pressure in the baryonic hot gas slows down the accretion 
with respect to the pressure-less "converted" dark matter 
component (especially in low density regions where the cool- 
ing time of the gas is extremely long). We speculate that the 
resulting delay in formation within the hydro simulation in 
this case could explain the delayed accretion of hydro satel- 
lites shown in Fig. [7] Nevertheless further studies are needed 
to confirm our hypothesis. 
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3.1 Orbital Parameters 

The largely different time evolution of R for twin haloes 
shown in Fig. [6]suggests that the radial position at any fixed 
redshift is not well-suited for directly comparing the orbits 
of the "twin" satellites. The peri-centric and apo-centric dis- 
tances (see Section [2,2[) . as well as the length of the semi- 
major axis of the orbit R avg = {Rperi + R ap o) /2.0, seem a 
better choice. 

Figure[T0]shows the ratio of the average distance (Ravg) 
for hydro and pure DM satellites as a function of the av- 
erage distance in the hydro run. In the pure DM case 
Ravg = Ravg/Rvir tends to exceed the corresponding value 
in the hydrodynamical simulations for satellites that are 
close to the center. In the outer region of the halo the trend is 
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Figure 10. C omparison of average distance of satellite "twins" at redshift 
z = 0. Each couple is color coded according to its DM observation time (i.e. 
the time at which they arc closer than 2 X R v i r from the main halo center). 



reversed and satellites have a smaller average distance in the 
pure DM case. This change in the radial behaviour seems to 
happen around R aV g > 0.85. We will use this (somehow ar- 
bitrary) threshold to separate the inner and outer behaviour 
of our twin subhaloes in the next figures. 

In this same plot, twin couples are color coded according 
to their observation time (i.e. the time at which they are 
closer than 2 x R V i r from the main halo center). We find 
a very interesting correlation between the ratio of R aV g in 
the DM and hydro runs and this observation time. Satellites 
with low t oos tend to live in the outskirts of the halo and are 
on average closer to the halo center in the DM run than in 
the hydro case. In contrast, "old" satellites (i.e. large t b s ) 
tend to live in the inner region of the halo and are closer to 
the halo center in the hydro run than in the pure DM. 

This result can be understood in the following way: in 
the outer region of the halo, dynamical friction is less im- 
portant, and the distance from the center of a given satellite 
mainly depends on its accretion time. On the other hand, 
in the inner regions dynamical friction is stronger in the hy- 
dro simulation (due to the central stellar body), dragging 
satellites toward the center in a more effective way com- 
pared to the pure DM case. This picture is confirmed by 
the comparison between the apo-centre distances for twins 
shown in Fig. 1111 Satellites in the outer region of the halo 
{Ravg,h y dro > 0.85) have a larger apo-center distance in the 
hydro simulation, while the situation is reversed for satellites 
further inside the halo. 



3.2 Mass loss 

The rate at which mass is removed from a given satel- 
lite depends on the balance between the density profile 
of the satellite and th e density profile of the central halo. 
iLibeskind et all l|20ld ) found that satellites in hydro simula- 
tions experience a lower mass los s due to their increas e d cen - 
tral density. On the other hand iRomano- Dfaz et all (|2010T ) 
found that the central region of hydro galaxies is depleted 
of satellites at a faster rate than in pure DM simulations, 
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Figure 11. Comparison of apo-centric distance of satellite "twins" at red- 
shift z = 0. The apo-centric in computed using a frozen potential modelled 
according to eq: 1. 



Figure 12. Scatter plot of relative retained mass of twins in pure DM 
and hydrodynamical simulation (combined results for G0,G1,G2) 



pointing to a faster and stronger mass loss in the hydro 
case. Figure [12] shows a scatter plot of the retained mass 
of all twins at z = 0. Satellites that lose less than 60% of 
their mass tend to retain more mass in the hydro simul ation, 
consistent with the findings of lLibeskind et all l|2010l ). But 
satellites that are heavily stripped tend to lose even more 
mass in the hydro simulation. There is also a clear corre- 
lation between the average position of the satellite and the 
amount of stripping. 

This result becomes clearer when comparing the mass 
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Figure 14. Same as figure fl3l but with twins color coded according to 
their mass at infall time in the hydro simulation (combined results for 
G0,G1,G2). 



loss with the orbital parameters of the satellite (R avg ), as in 
Fig. 1131 Satellites that live in the external part of the halo 
tend to be at larger distances from the center and retain 
more mass in the hydro simulation than in the pure DM one. 
On the other hand, satellites in the central region are closer 
to the center and more heavily stripped in the hydro sim- 
ulation, pointing to a more efficient dynamical friction and 
tidal stripping (in the central regions) i n the hydro simula- 
tion, i n agreement with the findings of iRomano-Dfaz et al.l 
(|2010h . 

In Fig. [14] satellites are color-coded according to their 
total mass (in the hydro simulation) at the time of infall. 
As expected from dynamical friction theory, more massive 
satellites are living at z = closer to center than low mass 
ones. This implies that satellites with large Mi n f a ii will be 
even closer to the center in hydro than in the DM simula- 
tions. 
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Together these results suggest a bimodal picture in 
which high mass satellite at the time of accretion are closer 
to the center and more heavily stripped in the hydro sim- 
ulation than in the pure DM, but the reverse for satellites 
with low Minfaii, which are on average further way from the 
center and less stripped in the hydro simulation. 



3.3 Danger to the galaxy disk 

In this section we assess if and how the different or- 
bital parameters and mass loss rates found for hydro mas- 
sive satellites influence results previously obtained using 
pure Nbody simulations (iFont et al.ll200ll; iRead et afll200Sl ; 
IVillalobos fc HelmillioOSl ; iKazantzidis et alj|2009h . 

Figure [15] shows a scatter plot of mass versus peri- 
centric distance for two different twin substructure popu- 
latio ns within host gal a xies G 0-G2. It is equivalent to Fig. 
1 of iKazantzidis et al.l (|2009l ) and we have used the same 
values for the scaling quantities Mdisk and Rd as in their 



work, namely M disk = 3.52 x 1O 1U M and R d = 2.82 kpc. 

The first substructure population in Fig.[l5]is comprised 
of all systems that have crossed an infall radius of r- ln { = 
50kpc from their host halo center since redshift z = 1. This 
selection is empirically fixed to identify orbiting satellites 
that approach the central regions of the host potential and 
are thus likely to have a significant dynamical impact on 
the disk structure. We assign masses to the satellites of this 
group at the simulation output time closest to the time of 
the first inward crossing of rinf, and then use the potential 
at this same simulation output to compute the orbit peri- 
center. Note that a single distinct object of this population 
may be recorded multiple times, as one subhalo may undergo 
several passes through the central regions of its host with 
different masses and peri-centers. Many of these satellites 
suffer substantial mass loss prior to z = 0. 

The second subhalo population consists of all surviving 
substructures at z = 0. The dotted line in Fig. 1151 encloses 
an area in the M su b — r pC ri plane corresponding to satel- 
lites more massive than 0.2Afdisk with peri-centers of r- pC ri < 
20kpc (r pcr i < 7Rd)- Subhaloes w ithin this area should be 
effective perturbers, and following Ka zantzidis et all (|2009l ) 
we refer to this area as the "danger zone" . We note, though, 
that we include this definition here for the purposes of illus- 
tration, only. 

By this criterion, Figure [15] suggests that satellites in 
hydro simulation are slightly less dangerous to disk stability 
than their pure DM counterparts. This can be interpreted on 
the basis of results shown in the previous sections. While the 
peri-centric distance does not substantially differ in hydro 
and DM runs, hydrodynamical satellites face a stronger mass 
loss. Since the danger for galactic disks is higher the more 
massive the satellite, this faster mass depletion for hydro 
satellites implies less dangerous perturbers for galactic disks. 



4 DISCUSSION AND CONCLUSIONS 

With this work we aim to provide a detailed description of 
the effects of baryonic physics on the properties of galactic 
satellites and their evolution with redshift. For this purpose 
we analyse three different cosmological simulations of galaxy 



10 k< 



"a 



DC 



0.1 



o o A 
° * 



DM z=0 
Hydro z=0 
DM 0<=z<=1 
Hydro 0<=z<=1 
danger zone 



0.01 



0.1 



1 



10 



M sat /M disk 



Figure 15. Scatter plot of peri-centric distance versus satellite mass. 
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50 kpc after 2 — 1 while unfilled symbols mark the properties of surviving 
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formation, each run twice: once as pure dark matter and an- 
other with the addition of gas physics, including gas cooling, 
star formation and feedback. For each run we create a com- 
prehensive catalogue of the subhalo population at z — 
which we trace back in time in order to study mass loss, 
dynamical friction and evolution of satellite orbital param- 
eters. Within this population we focus on a sub-sample of 
corresponding DM and hydro satellite pairs ("twins") and 
their individual evolutionary histories. 

Satellites are found to be more radially concentrated 
in the hydro simulat ion than in the pure DM one, confirm- 
ing earlier results bv lMaccio et al.l (|2006l ). This bias persists 
also for the twin population, even if slightly less pronounced. 
When we restrict our analysis to the twin sub-sample we 
find that hydro satellites tend to enter the virial radius of 
the parent halo later than the corresponding DM subhaloes, 
with an average delay of 0.7 Gyrs. This difference cannot 
be ascribed to a difference in the evolution of R V i r in the 
two simulations, and we speculate that it is instead related 
to the pressure support of the hot gas that acts against col- 
lapse, which, in low density regions, is not counterbalanced 
by cooling. Nevertheless further studies are need to confirm 
our hypothesis. 

Given the delay in the accretion time, the orbits of twin 
satellites are often weakly correlated. As a consequence, we 
find that the radial position at a given redshift is not suffi- 
cient to describe the satellite orbit. For this reason we define 
an average radius R avg for each satellite equal to the me- 
dian between the apocenter and pericenter distances of the 
satellite orbit, where these last two quantities are obtained 
by integrating each satellite orbit in a fixed potential resem- 
bling that of the halo. We find that the ratio of the average 
positions in the hydro and DM cases measured at redshift 
zero correlates with the accretion time of the satellite and 
its mass at that time. Moreover, we find that both the abso- 
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lute mass loss experienced by the satellite and the difference 
in mass loss in the hydro and DM simulations also correlate 
with the subhalo average distance at z = 0. 

We arrive at a final picture in which more massive satel- 
lites at the time of accretion are closer to the center and 
more heavily stripped in the hydro simulation than in the 
pure DM. The situation is reversed for satellites with low 
Minfaii that are on average further way from the center and 
less stripped in the hydro simulation. 

This bimodality can be understood in the following way: 
in the outer region of the halo, dynamical friction is less im- 
portant, and the distance from the center of a given satellite 
mainly depends on its accretion time. Since hydro satellites 
tend to be accreted later, they experience lower dynamical 
friction and mass stripping than their pure DM counter- 
parts. In the inner regions, on the other hand, dynamical 
friction is stronger in the hydro simulation (due to the cen- 
tral stellar body) and drags satellites toward the center in a 
more effective way compared to the pure DM case. Although 
the baryons provide a substantial glue to the subhaloes, the 
main halo exhibits the same trend. This assures a more ef- 
ficient tidal stripping of the hydro subhalo population, re- 
sulting in a larger mass loss. 

During the makin g of t his work, two other groups , i.e. 
iLibeskind et all l|201Ch and H omano- Dfaz et alJ ^OlOh . in- 
dependently pursued the same subject. Both studies com- 
pared hydro and pure DM simulations using an approach 
similar to ours, including a focus on corresponding pairs of 
satellites. The publication of Libeskind et al. mainly con- 
siders the radial distribution of the satellites and a com- 
parison of the retained masses. Like us, they find a differ- 
ence in the orbit s of tw ins. But in contrast to our approach, 
ILibeskind et~ai1 ||201Ch use the statistics of the radial or- 
bit position, rather then, e.g. determine the orbital parame- 
ters and the resulting average distance as we do. They find 
a larger mass loss for the satellite in the pure DM simu- 
lation and interpret this as a sign of the expected higher 
stability of the hydrodynamical satellite, which our results 
are only able to con f irm in the external region of the halo. 
iRomano-Di'az et al l (|2010h on the other hand, find an in- 
creased mass loss for the hydrodynamical satellite in the 
central region of the halo, in accordance with our results. 
They are also able to detect the final disruption of a satel- 
lite, and have shown that the life expectancy of the satellites 
in the hydrodynamical simulation is indeed shorter than in 
the pure DM simulation, as our findings suggest. 

We extend this result further in the last part of our 
study by investigating the possible impact of different orbital 
parameters and mass loss in hydrodynamical simulations, 
i.e. whether this translates into an increase, or a reduction, 
in the danger that these satellites pose for the stability of 
a possible stellar disk at the center of the parent halo (e.g. 
iKazantzidis et alj[2009t iMoster et aUl201Ch . While the peri- 
centric distance does not substantially differ in hydro and 
DM runs, hydrodynamical satellites in the central region 
face a stronger mass loss. As the danger for galactic disks is 
linked to the mass of the satellite, this faster mass depletion 
for hydro satellites leads to less dangerous perturbers for 
galactic disks, possibly easing the problem of the existence 
of thin disks in a Cold Dark Matter Universe. 

We emphasize, t hough , that, as correctly pointed out by 
iRomano-Diaz et al.l (|2010h , the effects of baryonic cooling in 



the center of dark matter (sub)haloes can be altered (if not 
reversed) for a more efficient feedback from stellar evolu- 
tion and possibly central super-massive black holes, which 
will expel baryons from the center and decrease the centra l 
concentration of the prime halo (e.g. |E1-Zant et ail 12004 
iMashchenko et al.ll2006l ; lGovernato et alj|20ld ). 

At this time, the direct comparison of our study of 
galactic satellites with observations is not possible. The 
Sloan Digital Sky Survey, which greatly contributed to the 
discovery and study of several Milky W ay satellites (e.g. 
IKoposov et"afll200Sl : Ide Jong et ai1l2010h , covers only a sin- 
gle patch of the sky. Future surveys like Pan-Starrs, which 
will provide a more comprehensive map of the (northern) 
sky, promise a better understanding of the properties and 
spatial distribution of satellite galaxies orbiting around the 
Milky Way. 
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